Tight bounds for the quantum discord 
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Quantum discord quantifies quantum correlations beyond entanglement and assumes nonzero 
values, which are notoriously hard to compute, for almost all quantum states. Here we provide 
computable tight bounds for the quantum discord for qubit-qudit states. In the case of two qubits 
our lower and upper bounds coincide for a 7-parameter family of filtered X-states, whose quantum 
discords can therefore be evaluated analytically. An application to the accessible information of 
the binary qubit channel is also presented. For the qubit-qudit state output by the circuit of 
deterministic computation with one qubit, nontrivial lower and upper bounds that respect the zero- 
discord conditions are obtained for its quantum discord. 
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The quantum discord [l], 0] has gradually become an- 
other important resource for quantum informational pro- 
cessing tasks besides the entanglement. Firstly, in cer- 
tain quantum mechanical tasks such as the deterministic 
quantum computation with one qubit (DQC1) 0, the 
quantum advantages can be gained Q with the pres- 
ence of quantum discord while the entanglement is ab- 
sent. Secondly the quantum discord is shown to be a 
better indicator of the quantum phase transition in cer- 
tain physical systems than the entanglement [j| . Thirdly, 
in addition to its original interpretation via Maxwell de- 
mon [g], the operational interpretations of quantum dis- 
cord via state merging Q establish firmly the status of 
the quantum discord as an essential resource amidst other 
concepts of quantum information. 

As a measure of the quantum correlation beyond en- 
tanglement, the quantum discord of a given state g of a 
composite system AB is [3, Q 

D A (g) ■= min y2piS(g B \i) + S(g A ) - S(g), (1) 

where S(g) = — Tr(g\og 2 g) denotes the von Neumann 
entropy and the minimum is taken over all the positive 
operator valued measures (POVMs) {Ef} on the subsys- 
tem A with pi = Tr(Ef-g) being the probability of the 
i-th outcome and Qbu — Tr A (Ef g)/pi being the con- 
ditional state of subsystem B. The minimum can also 
be taken over all the von Neumann measurements 
These two definitions are in general inequivalent even for 
qubits and our proposed bounds in Eq.© below apply 
to both of them. When the measurements are made on 
subsystem B the quantum discord Dsig) can be defined 
similarly and is in general different from D A (g). 

Quantum discord assumes nonnegative values and 
zero-discord states are relatively well understood: 
D A (g) = if and only if there exists a complete or- 
thonormal basis {\k}} for the subsystem A together with 
some density matrices g? for the subsystem B such that 



8 = SfcPfcl^)(^l ® J?fcL, Recently various methods to de- 
tect zero discord [8|, [9| have been proposed for a given 
state as well as for an unknown state Besides its 

initial motivation in pointer states in measurement prob- 
lem vanishing quantum discord is also found to be 
related to the complete positivity of a map II ll and the 
local broadcasting of quantum correlations |12j | . 

Unfortunately zero-discord states are of zero measure 
[l3j and the nonzero values of the quantum discord are 
notoriously difficult to compute because of the minimiza- 
tion over all the possible measurements. There are only a 
few analytical results including the Bell-diagonal states 
14 1 and rank-2 states 11511 . in addition to a thorough 
numerical calculation [16| in the case of von Neumann 
measurements. For two-qubit X-states, since there are 



counter examples Yjj, [18[ for the algorithm proposed in 
(ljjj . the evaluation of their quantum discords remains 
a nontrivial task. It is therefore desirable to have some 
computable bounds for the quantum discord so that the 
question of how large or small the quantum discord can 
possibly be, e.g., in the DQC1 circuit, can be answered 
more reasonably. 

In this Letter we shall provide computable tight 
bounds for the quantum discord of qubit-qudit states. 
For a family of two-qubit filtered X-states with 7 pa- 
rameters up to local unitary transformations (LUTs) our 
lower and upper bounds coincide so that their quantum 
discords can be evaluated analytically. Also we present 
an application to the accessible information of the binary 
qubit channel for which our lower and upper bounds can 
coincide. For the quantum discord of the qubit-qudit 
state output by the DQC1 circuit we derive nontrivial 
lower and upper bounds, which qualitatively agree with 
the zero discord conditions [1, [lO] comparing with the 
estimation in 

To present our main result we need some notations. 
First of all, we denote by g a qubit-qudit state with the 
reduced density matrices g A and gs for the qubit A and 
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the qudit B respectively. Let x be the Bloch vector for 



qa with its length given by x 



2Trg 2 A - 



1. Without loss 



of generality we assume that the reduced density matrix 
qa of qubit A, on which the measurement is performed, is 
invertible since otherwise we have a product state, which 
has zero discord. Therefore the following filtered density 
matrix is well defined: 



1 



1 



(2) 



Secondly, let ay, (// = 0,1,2, 3) be the identity matrix 
and three standard Pauli matrices. We associate with 
every qubit-qudit state g (or g) a positive semi-definite 
two-qubit operator 



2Tr BlB2 [(l-V£)g®g] 



(3) 



in which Vff 



= Y^ij stands for the two-qudit 

swapping operator. Then the equation det(Q e — qrj) = 0, 
where r\ = diag(l, — L— L— 1), has four real solutions 
Qi > Q2 > 93 > 94 120, \2M- It should be noted that 
each qt is as readily computable as the concurrence of 
p, which equals to max{0, ^fq^ + y/qj ± yfql — Jq~y\j2 
22|, in the case of two-qubit states for which we have 
Q e = RgT/R^ with R e being the 4x4 correlation matrix 
defined by [R 8 ]^ v = {^n®^f)g- These values qi(g), show- 
ing explicitly the dependence on g, are invariant under 
an arbitrary Lorentz transformation (LT) L, satisfying 
LrjL T = r], that brings Q g to LQ e L T . The local filter 
\/l — x 2 / \/2g~A acting on qubit A brings Q e to (1 — x 2 )Qg 
and induces an LT to Q g , from which it follows that 
(1 - x 2 )q 2 (g) = q 2 {g). 

Thirdly, we denote by Q| x3 the 3x3 matrix obtained 
from the 4x4 matrix — Qg by deleting its first row and 
column. Let ti(g) be the largest eigenvalue of Q| x3 cor- 
responding to the eigenvector e. We define fhg to be the 
unit vector along the direction yl — are± + eji , where 
ej_ = e — e]| and ej| = x(x ■ e)/x 2 . In the case of x = 
we define fhg 



e. By measuring the observable frig ■ a 
on the qubit A we obtain a suboptimal value DA(g\rrig) 
for the quantum discord as in Eq.(JT|) without minimiza- 
tion. For a two-qubit state g, since TtbQ = ctq/2, we have 



)? x3 
'S 



TgT~ where Tg is the 3x3 correlation matrix 



for g defined by [Tg] ab = (a a (g) a b )g (a, 6 = 1,2, 3). 

Lastly, we need to introduce an increasing convex func- 
tion co(l — z 2 ) of z where 



CO(Z) = < , 2 

^ l0 S2 TTI: 



if z > 0, 
if z < 0, 



(4) 



and = log 2 _ 1^5 i og2 1^£. For any 

density matrix g it holds S(g) > co(2Trg 2 — 1) since we 
always have S(g) > — log 2 Tig 2 and if Tig 2 > 1/2 we 
have S(g) > /i(2Trg 2 — 1) which can be read off from the 
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FIG. 1: (Color online) Comparing the upper and lower bounds 
with the quantum discord minimized over von Neumann mea- 
surements for 10 randomly chosen 2-qubit states of rank 4. 



information diagram between the entropy and index of 
coincidence [23j |. 

Theorem For a qubit-qudit state g, with subsystem A 
being a qubit, the following bounds hold 

co(C) + S(g A ) - S(g) < D A (g) < D A (g\A s ), (5) 

where C = 2Tvg 2 B — 1 + q 2 (g). In the case of two qubits, 
the lower and upper bounds coincide if q 2 (g) — t\(g). 

A detailed proof can be found in 24] and an out- 
line is given in what follows. The upper bound holds 
true by definition. The key to p rove the lower bound is 
the Koashi- Winter relation [25(: .Da ( Pas) + S(qab) — 
S(g A ) = Ef(p_rc), where Ef(qbc) is the entanglement 
of formation [26( of the dx2d state gsc = ^A^i^ABC 
with \ip) abc being a purification of gAB in a 2 x d x 2d 
system ABC. The lower bound follows from Ef{qbc) > 
co(l — C BC ), where Cbc is the concurrence of gBC, and 
an evaluation of the concurrence C 2 BC = 1 — C by apply- 
ing Theorem . 
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Furthermore the ensemble of qbc 
induced by measuring fhg ■ a on qubit A is optimal for the 
tanglcTsc' 01 Qbc whose value 2(1 — Trpp) — (1 — x 2 )U (p) 
can be readily obtained by applying Theorem 1 in [27| . 
Since h(l — z)isa concave function of z we have bound 
D A (g\rng) < h(l - t B c) + S(g A ) - S(g AB ), which co- 
incides with the lower bound if q 2 (g) = ti(g), i.e., 
tbc — Cbc> since C > in the case of two qubits. 

For two-qubit states we have carried out some numer- 
ical calculations and a comparison with our bounds is 
shown in Fig.l, where for convenience the upper bound 
is taken to be h(l — tbc) + S(g A ) ~ S{gAB), which is 
slightly weaker than Eq.©. For about 70% of 10 5 ran- 
domly chosen states of rank 4 the differences between our 
bounds and the values of quantum discord obtained by 
minimization over von Neumann measurements lie in the 
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range of ±0.01 with a maximal difference about ±0.1. We 
note that for some states our lower bounds are negative 
and therefore are trivial. On the other hand there exist 
several nontrivial families of states, for which the upper 
and lower bounds coincide so that analytical expressions 
of quantum discord for these states can be obtained. In 
these cases the two definitions of quantum discord be- 
come identical. 

Example 1: Bell-diagonal states fl^i l. A Bell-diagonal 
state has a density matrix qbs — \ C/i ^® ^ with 
being real and |c M | < Co = 1. It is clear that qbs = Qbs 
and Q eBS = diag(l, — cf , — c|, — c|). As a result we have 
o_2{qbs) = ti(gBs) = max{c 2 , c§, c|} so that upper and 
lower bounds as in Eq.([5]) coincide. 

Example 2: Rank-2 states )!&] ■ Every two-qubit state 
of rank 2 admits a 3-qubit purification \i)))abc with the 
reduced density matrix qbc being a two-qubit state of 
rank 2. As a result r^c = C BC so that the upper and 
lower bounds as in Eq.([5]) coincide. 

Example 3: X -states. In many scenarios A-state arises 
as the two-particle reduced density matrix as long as 
there is a certain symmetry of the physical system. In 
general a 2-qubit X-state 



Qx 



( Qoo O3 ^ 

011 012 

021 Q22 

\ 030 033 / 



(6) 



has 7 independent real parameters. Since the quan- 
tum discord is invariant under LUTs, we can assume 
without loss of generality that 003 and 012 are real 
and therefore we have in fact only 5 real parameters, 
which can be conveniently taken as those nonzero en- 
tries of the correlation matrix R ex , i.e., x = (03 ® o~o) ex , 
V = (0-0 ® cr 3 ) gx , and s a = (<r a <g> a a ) gx (a = 1,2,3). 
Since Q e = Rg-qR^ holds for qubits it is easy to see that 
four solutions to the equation det(Q 0Y — qw) — are s 2 2 
and 4(^/000033 ± \jQwQ22)' 1 ■ Furthermore the 3x3 corre- 



lation matrix Tg x is diagonal with entries Si^/vl — x 2 
and (S3 — xy)/(l — x 2 ). Thus as long as 



I V 000 033 - y/Qll&Vl] < |0O3| + 1 012 1 



we have qziQx) = H{Qx) = max{sf 2 }/(l " x2 ) so that 
the upper and lower bounds of the quantum discord co- 
incide. We denote by Qx a the X-state that satisfies the 
condition Eq.© and the optimal observable to measure 
is (Ti if s\ > S2 and o~i otherwise. Explicitly, 



Da(qx c ) = h(y 2 + max{s? j2 }) + h{x 2 )- 
1 ± s 3 , f (x ± y) 2 + (si T s 2 ) 



± 



(l±s 3 ) 2 



(8) 



Example 4-' Filtered X -states. Take an arbitrary 2- 
qubit state g satisfying i) 92(0) = ^1(0), i-e., the up- 
per and lower bounds for Da{q) coincide, and ii) qa is 



proportional to identity, for example the state gx c with 
x = as in example 3, and take an arbitrary invertible 
Hermitian operator Fa acting only on qubit A. Then for 
the filtered state Qp cx FaqF\ we have g = gp from the 
definition Eq.([2]) so that gF has also a coincided lower 
and upper bound. Because cFa and Fa represent the 
same filter for any nonzero real c, there are 3 indepen- 
dent real parameters for Fa- Thus we have a family of 
filtered X-states with 7 parameters for which the quan- 
tum discords can be evaluated analytically. 

Example 5: Accessible information. A quantum com- 
munication channel {pk,gk} is defined by the action of 
sending states {0fc} with probabilities {pk} to a receiver 
who can perform any possible POVM to gain the infor- 
mation of {pk}- The accessible information I acc is the 
maximum of the mutual information of the joint prob- 
ability distribution {p a k = PkREaQk} over all POVMs 
{E a }. If we introduce a bipartite state gAB — SfcPfePfe® 
|fc)(fc| where {|fc)} is a complete orthonormal basis then 
we have D b {qab) = and Da(qab) = x{qa) - Iacc 
where x(qa) '■— S(qa) — J2PkS(gk) is the Holevo bound. 
For qubit channels the theorem provides an upper bound 
Iacc £ S(qb) — co(C) as well as a lower bound given by 
the measurement along rh g for the accessible information. 

Consider a binary qubit channel {Pk, 0k} 2 = i and de- 
note by a and b the Bloch vectors of qubit states 01 and 
02, respectively, and c± = p\a±p2b. The lower bound 
for I acc provided by rh e coincides with the lower bound 
(M(i)) in 
qubit state g - 



The 4x4 correlation matrix of the two- 
=1 PkQk ® \k){k\ reads 



El 



Rg = 



1 5 







(9) 



where S = p\ — p2- We note that Q e — R g j]R^ is of rank 
2 and therefore four solutions to the equation det(Q g — 
qrj) = are {0, 0, (1 — 6 2 )A±}, where 



A, 



l-a-b± v /(l-a 2 )(l-6 2 ) 



(10) 



from which we obtain 52(0) = (1 — <5 2 )A_. The 3x3 
correlation matrix Tg of g is of rank 1 and the largest 



(7) eigenvalue of TgTj is t^g) = (1 - S 2 ) 2 X + X 



7(i 



r.2 \2 



By solving 52(0) = ^1(0) we obtain p\{l — a 2 ) — p\{\ — 
b 2 ) in which case the upper and lower bounds for I acc 
coincide, meaning that the optimal POVM for I acc is a 
von Neumann measurement (along the direction m g oc 
c_). This proves explicitly, for the first time as far as 
we know, that the optimal von Neumann measurement 
given in [28( is also optimal among POVMs. 

Example 6: DQC1 states. Let us consider the qubit- 
qudit state output by the DQC1 circuit, aiming at calcu- 
lating efficiently the trace of a unitary operation U acting 
on the qudit B, with a density matrix 3| 



Qu = TTi 



1 

27/ 



aU 



aW 
I 



(11) 
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where d = 2" and < a < 1. Let {e^*, |a)}f =1 be the 
eigensystem of U and obviously g u = h Tla^a ® l a )( a l 
is separable where <j) a = (1 + cm a • c?)/2 with n a = 
(cos</> a , sin </> a , 0). The nonvanishing quantum discord of 
g u , which has been estimated in [4j, is argued to be re- 
sponsible for the quantum speedup. Recently it is found 
that [1, [l^ D A (g u ) = if and only if a = or /3 = 1 
where /3 = (d + \TrU 2 \)/(2d). The estimation in 
though respects the trivial condition a = 0, is insen- 
sitive to the condition /3 = 1, i.e., a nonzero value of 
the quantum discord is estimated in this case, in which 
U 2 oc I and U can be in the same time typical, i.e., 
til := |TrJ7|/d ~ 0. Our Theorem provides nontrivial 
bounds that respect both two zero-discord conditions. 

For simplicity we shall assume u% = first, i.e., E n a — 
0. In this case we have g u = g u and S(TibQu) = 1- It 
turns out that Q Bu — diag{2(l — 1/d),— Q 3 g * 3 } where 

QlT = IF Eti n a fi a has eigenvalues ^{0, /?, 1 - /3}. 
As a result we have qi{Qu) = 2a 2 (3/d (/3 > 1/2) which 
is also the largest eigenvalue of Q 3 e * 3 with eigenvector 
m = (cos (f>, sin 0, 0) where e l2<#> = Tri7 2 /|Tr[/ 2 |. From 
Eq.© and £ = 2(1 + a 2 /3)/d - 1 < for d > 4 it follows 

log 2 — ^-/i(« 2 ) < A^) < h (a 2 p)-h(a 2 ). (12) 
1 + arp 

In deriving the upper bound we have used D + 

M" 2 ) < iEa^a 2 ^;^) 2 ) < ^E„(m-n Q ) 2 ) in 
which the first inequality is valid for arbitrary u% and 
becomes an equality in the case of Ui = 0. Thus the 
upper bound for Da(q u ) given in Eq. ([T2"T) holds for all 
unitary U and reaches its maximum h(l/2) w 0.6 at a — 
1 and fj = 1/2. It is obvious that in the case of a = 
or in the case of pure qubit a = 1 and j3 = 1 the lower 
and upper bounds coincide so that the quantum discord 
vanishes. 

To conclude, in view of the hardship of computing the 
nonzero values of quantum discord, the computable tight 
bounds provided here for qubit-qudit states should be 
useful in further quantitative studies of the relation be- 
tween quantum discord and phase transitions, quantum 
speedups, and so on. Notably our bounds enable us to 
evaluate analytically the quantum discords of a family 
of filtered X-states with 7 parameters up to LUTs (for 
comparison there are 9 parameters for a general 2-qubit 
state) and to estimate the quantum discord in the DQC1 
circuit more reasonably. Also our bounds are applicable 
for the classical correlation and accessible information. 
Though we have restricted to qubit-qudit systems, for 
which the concurrence Cbc can be evaluated exactly, the 
bounds for the quantum discord Da{qab) of a general 
bipartite state, which might not be specially tight, can 
be obtained in a similar way via the bounds for the en- 
tanglement of formation Ef(qbc)- 
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SUPPLEMENTARY MATERIAL: PROOF OF 
THE THEOREM 

Let be the eigenstates of a given qubit-qudit 

state qab with eigenvalues = (ipilipi)- Then \iP)abc = 
Si=i IV'i) as ® \i)c is a purification of qab in a 2 x d x 2d 
system ABC. Then g BC = Tr^l^^lxBc is of rank 2 
and is supported on the 2-dimensional subspace spanned 
by an orthornormal basis 

W)bc=a (k\-^\ip) A BC, (A: = 0,1). (13) 

y/QA 

Any state xbc supported on the subspace spanned by 
{|0fc)sc} can be expanded xbc = \ S M r n a ^ C wit b the 
help of four generalized Pauli operators 



(14) 



k,k'=0 



where coefficients r M (r = 1) are real (fi = 0,1,2,3). 
For example we have qbc = {°~o C + % ■ <j bc )/2. From 



the fact Tr^cr 



bc _ 



2Ts:a{p'hQab) h follows that xb 



?> Tr a (c^ab) and 



2(1 - Tr X |) = 2Tr[(l - V&)xb ® Xb] 

3 

= 2 ^ Tr[(l - V5f )aj ® cr^gAB ® Sab]^ 

LLM — 

3 



(15) 



in which Qg has been defined in Eq.©. The proof of the 
theorem is divided into the following four steps. 

Lemma i) Ef(qbc) > co (l — Cbc)> u ) ^bc = 1 ~ 
iii) t bc = 2(1 - Trg|) - (1 - x 2 )h(g), and iv) For a 
two-qubit state g^s the bound 

D A {QAB\mg) < h(l - T B c) + S(qa) - S(qab) 

holds true. The entanglement of formation, concurrence, 
and tangle of qbc are defined to be, respectively, 



E f {qbc) = mm VpiSVf), 



(16) 



C BC = min V K J2(1-Tr(7rf)2), (17) 



T~BC 



mm 

{p.,^f c } 



\ " 



2 K (1-Tr(7r 



B\2\ 



(18) 



where the minimization is taken over all possible en- 
sembles qbc = ^iPi^f ' wrt h 7r l SC being pure and 



BC • 



irf = Tr C 7r. 

Proof, i) Let {pi,n B } be the optimal ensemble for 
Ef(qbc) and we obtain 



Ef{qbc) 



PiS(*f) 



> 5> CO (2Tr(^f) 2 - 1) 



> co(l-Ci c ), 



(19) 



in which the first inequality is due to S(g) > co(2Trg 2 — 
1) the second and third inequalities are due to the fact 
that co(l — z 2 ) is a convex and increasing function of z, 
respectively. 

ii) Since qbc is of rank 2, its concurrence Cg C — 
2(1 — Trg^) — (1 — x 2 )q2((>) can be exactly evaluated ac- 
cording to Theorem 4 in 21(, where 92(g) is the second 
largest solution to the equation det(Qg — q-q) = 0. Since 
the local filter yl — x 2 / ' \J2qa acting on qubit A brings 
Q e to (1 — x 2 )Qg and induces an LT to Q g we have 
(1 — x 2 )q2(g) = 52(g) and consequently Cfj C = 1 — L. 

iii) For an arbitrary state xbc = ^J2 fJ , r ^ (J ^ C 1 with 
Xb and xc denoting the reduced density matrices, we 
introduce a 4 x 4 matrix M whose elements are given by 

1 - Trxi - Tr Xc + T^Xbc 
= Tr[(l - Vg )(1 - Vg) X BC ® Xbc] 

= 1 £ Tr[(l-^)(l-^)< c 



1 <^ c ]r M rv 



' /n,i/=0 
3 

E M ^ 

fj,,v— 



v> 



(20) 



According to Theorem 1 in [27J we have r^c = 1 — Trg^ 
Trg^+Trgg C ;-|-(l— rc 2 )f where i is the smallest eigenvalue 
of the 3x3 matrix M 3x3 obtained by deleting the first 
row and column of the 4x4 matrix M defined above. It 
turns out that M = Qg + |(1 — a)rj where 



a = Tr[(Vif - Vg)\<t> )(cf> \ ® |0i)(0i|] = 



Trg| - Trg 2 , 
1-Trg| c ■ 



Thus t = — ti(g) + (1 — a)/2 which proves Lemma iii. 

iv) Let {p±,n BC = {a BC + e± • ct bc )/2} be an op- 
timal decomposition of qbc for the tangle tbc and we 
have e+ — e_ cx e, the eigenvector of TgTj correspond- 
ing to its largest eigenvalue ii(g), and x = p+e + + 
p_e_. If trig is the unit vector along the direction 
yl - x 2 ej_ + eli then we have p± = (1 ± x ■ rfig)/2 and 
p±e± = Tr(ayfgA(l ± m§ ■ <?)\/&a) /2- Therefore if wc 
measure the observable mg • ct on the qubit A we obtain 
outcome ± with probability p± leaving qubit B in the 
state q B \± = Tta(1 ±£- mg)gAB/(2p±). Thus 



= Ep±^±) 

= 5> ± M2(Tr(7rf) 2 -l) 

< ft(l-TB C ) 



(21) 



in which the second equality is due to the fact that 7r^ 
is a qubit state and the inequality is due to the fact that 
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h(l — z) is a concave function of z. It follows immediately 
that 

D A (g\m~ e ) = ^2p±S(it±) + S(g A ) - S(qab) 

< h(l-TBc) + S(g A )-S(g AB ). (22) 

Proof of the theorem. The upper bound is triv- 
ially true. Because of the Kaishi- Winter relation [25| 
D a (q A b) + S(g AB ) - S(g A ) = E F (g B c), Lemmas i 
and ii lead to the lower bound in the first statement. 



In the case of two qubits, since C BC < 1 leading to 
C > 0, the lower bound of D A (g) in Eq.© becomes 
h(C) + S(g A ) — S(g A B)- From Lemma iv we see that 
h(l — tbc) + S(gA) — S(g A B) is an upper bound of 
D A (g\rfig) and therefore also an upper bound of D A (g). 
If 12(g) = t\(g) then tbc = C BC = 1 — C because of 
Lemma iii, which means that the upper and lower bound 
coincide. 



